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Abstract 

' The non-conserved t/)^ model defined by a Langevin equation with external 

non-white noise is studied by means of the Dynamic Renormalization Group. 
The correlation time of the noise changes the critical point location but does 
not affect the critical exponents up to order e^. The same effect is obtained 
when the correlation length of the noise is considered. These results are shown 
to be in agreement with previous numerical simulations. 
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I. INTRODUCTION 



Langevin equations in spatially-extended systems have been extensively used in the stud- 
ies of equilibrium and nonequilibrium phenomena where fluctuations are important. These 
equations are stochastic partial differential equations in which the fluctuations are introduced 
through a noise term, whose stochastic properties are chosen according to the physical sit- 
uation. Fluctuations of internal (thermal) origin represent microscopic degrees of freedom, 
and since they evolve in spatial and temporal scales much shorter than those of the gross 
variables of the system, they are assumed to be uncorrelated (delta correlated) in space and 
time. In this case they are modelled by a gaussian-and-white noise process in space and 
time. On the other hand, if the origin of the fluctuations is external to the system, the 
possibility of a noise with some spatial or temporal structure has to be considered. 

Langevin equations have been used in many different scenarios, one of which has been 
for many years the study of Dynamical Critical Phenomena The methodology used 

there is a suitable extension of the Renormalization Group techniques developed for the 
study of static critical phenomena The technique has also been generalized to study 

nonequilibrium systems presenting some kind of scale invariance or criticality 

In this paper we will follow this methodology as it was implemented in Refs. [^,0 to 
study the nonequilibrium critical properties of a Langevin-equation model with a non-white 
(colored) noise. 

Among the different models that can be considered we have chosen the simplest one, 
in which the order parameter is not conserved. This model is useful for describing order- 
disorder transitions. In this non-conserved case the equation of motion of the fleld variable 
is 



where x is the position in a d-dimensional space, F[%1)\ is the Ginzburg-Landau free-energy 
functional 



gV'(x,t) 
dt 



5F[%Ij] 



+ e(x,t) 
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F[i;] = I d'x f^Ki;' + ^/x I 1^ • (1-2) 

and the correlation of the noise is 

(e(x, t) e(x', t')) = 2D 5(x - x') 6{t - t'), (1.3) 

where D = ksT to ensure the correct equihbrium steady state. This is known as model A 
in the notation of Ref. 0. This model comes from a coarse-graining procedure apphed to a 
spin model 0. Its behavior can be easily understood by means of a mean-field analysis of 
the potential given by the Ginzburg-Landau free energy. When r is positive this potential 
has a minimum at zero field, which corresponds to a coarse-grained field equal to (i.e. 
the spins are randomly up or down). A negative value of r leads to a potential with two 
nonzero symmetrical minima (the spins have a preferred direction, either up or down). The 
role of noise in this scheme is a disordering one: when r is negative and the intensity of the 
noise increases, the potential barrier between the two minima is more easily surpassed by 
the spins. Thus, although the potential shape does not change, the coarse-grained field will 
eventually become zero (for a high enough value of the intensity of the noise): the spins will 
be disordered because their thermal energy (i.e. the noise intensity) is so high that they do 
not see the potential barrier. 



Our aim in this paper is to study the model (|1.1[ - |1.2|) in the special situation characterized 
by the fact that the noise term .^(x, t) is non-delta correlated. In this case there is no 
fluctuation-dissipation relation and the steady state is not the equilibrium one. Hence we 
are facing with a pure non-equilibrium problem. In a first approach to this problem we will 
assume that the correlation is of the Ornstein-Uhlenbeck type: 

< e(x, t)e(x', t') >= - e--^6\^ - x') (1.4) 

r 

This is a suitable way of modelling the noise in the case when it is external, r is its correlation 
time. It can be seen that in the limit r — > the noise becomes white and the equilibrium 
situation is recovered provided D = ksT. Naively, one might say that an increasing value 
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of T makes the noise "softer", taking it farther from the white character (given by r — > 0) 
and reducing its effective intensity. This effect has indeed been found numerically p] and 
approximate techniques borrowed from stochastic theory have been developed in order to 
understand this phenomenon Q]. In this paper we perform dynamic Renormalization-Group 
(DRG) calculations on this model in order to study the effects of the correlation time of the 
noise in the non-equilibrium critical behavior of the system. 

This paper is organized as follows: in next Section we present a preliminary scaling anal- 
ysis to obtain some of the critical characteristics of our model. In Sect. |T| the Dynamic 



Renormalization Group results are presented. In Sect. |V| we summarize the results corre- 
sponding to a noise colored in both space and time. The comparison between our theoretical 
results and those obtained from simulations is given in Sect. 0. Technical details of the 
DRG calculations are presented in the Appendix. 



II. PRELIMINARY ANALYSIS 

A glimpse of the relevance of the different terms in ( |1.2| ) at the critical point can be 
obtained by studying the effect of a scale transformation on the dynamical equation ( p..lD . 
Let this scale transformation be defined by 

x = bx; t = Vi- ^ = 6"'^ ^, (2.1) 

with 6 > 1, which means that one is focusing on long distances (characteristic of the critical 
point), a is related to the equilibrium exponent rj via a = |((i — 2 + ?7)|0]. Introducing these 
changes into Eq. ( |1.1| ) and making use of the expression of the Ginzburg-Landau functional 
( |1.2|) one can see that the dynamical equation becomes 

^ = -rh' ^ + 11 - u V-^" + 6"+" e (2.2) 

dt 

Now we define a rescaled noise C, = b"'~^^^, which has a new intensity 

D = D 62a+^-'^ (2.3) 
4 



and a correlation time r = t/V. Moreover, in order to keep the shape of the dynamical 
equation after the scale transformation, we also define the following rescaled parameters: 

f = rh\ fi = fi b'-\ u = u h'-"^" (2.4) 

Now we place ourselves in the critical region of the system (which in a first approximation 
is supposed to be located near r = 0) where, in the absence of the nonlinear term (u = 0), 
the equation remains invariant under the scale transformation if the exponents a and z are 
chosen to be 

d-2 

a = ao = {r]o = 0) 

z = zo = 2 (2.5) 

These are the trivial exponents of model (|1 . ■ 2|) with colored noise, which coincide with 
the ones corresponding to the white-noise case 0. They are valid when the nonlinearity of 
the model does not play a role in the critical behavior of the system, and this occurs beyond 
a critical dimension which we can now determine. As we have seen in ( |2.4|) , the nonlinear 
coefficient u transforms upon rescaling with an exponent z — 2a. By making use of the trivial 
values of z and a, one finds that this exponent happens to be zq — 2ao = 4 — d = £. Thus 
when c? > c/c = 4, the nonlinearity decreases under rescaling (6 > 1), so that u is irrelevant. 
On the other hand, for d < dc, u grows under rescaling, and hence it is now relevant. This 
analysis is independent of the character of the noise, i.e. it does not depend on whether 
the noise is white or colored. In fact, since z > and f = r b~^, the correlation time r is 
always an irrelevant parameter, so that one should not expect any effect of it in the universal 
properties of the model. However it has, as we will see, infiuence on the system. 



III. DRG RESULTS 

The DRG is going to be implemented in Fourier space, so the first thing we must do is 
to transform our Langevin equation ( |1.1[ ) conveniently. Let us define the Fourier transform 
of a field variable ip as: 



(3.1) 



where k is a vector in a d- dimensional (Fourier) space. According to these definitions, the 
Fourier transform of Eq. ( |1.1| ) can be written as 

-uGo{k,iu) / V(ki,u;i) ■?/'(k2,a;2) ■?/'(k - ki - k2,a; - co-i - 0^2) (3.2) 

Jki uJi J k2 UJ2 

where /^^ stands for ^^^^^+1 / d'^k du. is the zeroth-order aproximation (in u) to the field 
in momentum space: 



^°(k,o.) = Go(A;,u;)e(k,u;) 



(3.4) 



and Go is the propagator: 



Go{k,u) 



(3.5) 



r + /iA;^ — iuj 

On the other hand, by using definition ( pj.l| ) in Eq. (O) we can calculate the correlation of 
the noise in Fourier space: 



<e(k,^)e(k',^') >= (2 



TT 



1 + cuV^ 



5'^(k + k') 5{uj + cu') 



(3.6) 



- kg 



Diagrammatic version of Eq. ^S.Sj ) 
Fig. 1 

DRG will now be applied to Eq. ( p.2| ). In order to simplify this procedure, a diagram- 
matic notation is very convenient. The diagram version of Eq. ( |3.2| ) is shown in Fig. 1. 
Thick lines stand for the field, whereas a thin line represents the zeroth order aproximation 
of the field and a vertex stands for the two integrals (over ki and k2) and u. The propagator 
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is represented by a thin line with an arrow. As can be seen in this figure, momentum through 
a vertex must be conserved. 

In this work we have decided to carry out a DRG analysis using the standard momentum 
shell integration scheme. While details of this method are well documented in Ref. [0, here 
we briefiy review the procedure, which consists of the two following steps: 

a. Momentum Shell Integration. The first step of the DRG is to eliminate the short 
wavelength modes (i.e. integrate out modes in the momentum shell Ae^' < |A;| < A, where 
/ > 1 and A is the momemtum upper cutoff), since we are interested in the long wavelength 
behavior of the system. 

b. Space Rescaling. After the momentum shell intergration step, we rescale space in 
such a way that the full momentum space is recovered. This leads to the differential equations 
satisfied by the running coupling constants. A fixed point analysis of these equations will 
then reveal the scaling properties of the long wavelength correlations which we are seeking. 

Explicit details of the DRG calculations for our model are summarised in the Appendix. 
The DRG results are obtained from the analysis of the differential flow equations 
defined as the infinitesimal variation of the parameters when the renormalization step / is 
very small. From ( [A 61 ), they are simply: 

dr 

— = zf + 3 u KJl — jj,T) (3.7a) 
dl 

dzL 

— = {z- 2a)u - 9uu K4 (3.7b) 
^ = - 2) /i (3.7c) 
^ = {2a + z-d) D (3.7d) 

where the following parameters have been defined: 

r uD 
f=-, 3.8 

The fixed point of this transformation is the one at which these derivatives are zero. 
Physically it represents the critical point of the system. Therefore, by studying its position 
as a function of r we will be able to study the infiuence of the time correlation of the noise 



in the disordering transition induced by the noise intensity. By imposing this invariance on 
Eqs. ( p.7c| ) and ( p.7d| ) one finds that the values of the exponents a and z up to 0{e) are 



the trivial ones already obtained in Sec. 0. Invariance of the static parameters (Eqs. ( p. Taj) 
and ( |3.7b| )) shows that the fixed point is given by: 

r = -|(l-.-r). r.- = ^ (3.9) 

For r = (white noise) the critical value of r is negative, which is reasonable: in the presence 
of an additive (disordering) noise, the system will be disordered even for a small negative 
value of r. When r 7^ the critical value of r is nearer 0, which means that the disordering 
effect of the noise will somehow become diminished because of its correlation in time. Thus 
the effect of the time correlation of the noise is a "softening" one, as expected and explained 
above. It is worth noting that in the white-noise limit our results coincide with the known 
values appearing in the literature 

Concerning the evaluation of the critical exponents of the system, this can be done in a 
simple way once the recursion relations are known. By making use of standard techniques 
involving the evaluation of the eigenvalues of the recursion-relation matrix it is easily 
seen that there is no contribution of r (up to first order) to any of the exponents. Hence one 
can conclude from this analysis that a correlation in time of the fiuctuations affecting a 0^ 
non-conserved model infiuences the position of the critical point of the system (as observed 
numerically 0), but not its universal properties (through its critical exponents). 



IV. COLORED NOISE IN SPACE AND TIME 

A more realistic assumption in relation to external noise is the existence of a non-delta 
correlation also in space. Intuition tells us that the infiuence of a non-zero correlation 
length of the noise will be qualitatively the same as the effect of the correlation time, i.e. 
an ordering one. This was checked numerically and also explained theoretically in Ref. 0. 
A DRG argument concerning this behavior can be made in a straightforward way following 
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the procedure previously described. The noise is chosen to be a generahzation of Eq. ( [L.4|) 
with no delta correlation in space. This can be done in a simple way by assuming that it 
obeys the Langevin equation 

t) = - i (1 - A e + - Vir, t) , (4.1) 

T T 

where ri{f,t) is a gaussian white noise with correlation (|1.3| ). This is an extension of the 
Ornstein-Uhlenbeck process, where the laplacian term takes into account the coupling of the 
field at different points [|lOl , so that A is the correlation length. The correlation of this noise 



in Fourier space can be seen to be 

< e(k, u') >= {2nY^'—-——^-^ + k') 6{uj + u'), (4.2) 

(1 + X^k'^) + cu^r^ 

which can be compared to Eq. ( p.6|) . 

The scaling analysis of Sec. |I| applied to this new situation shows that A changes as 
A = A/6, so that it is an irrelevant parameter (as r). Nevertheless we aim to analyze, as in 
the case of the correlation time, its nonuniversal effects on the system. 

The contribution of A to the diagrams which renormalize r and u can be easily obtained. 
The new results for the corresponding differential flow equations are 

^ = zf + 3 u Ki ^ (l (4.3a) 

^ = (z- 2a)u - 9uu ^ ^ (4.3b) 

dl ^ ' ' (1 + A2)2 ^ ' 



which should be compared to Eqs. ( |3.7a| - pT5|) . These new contributions of A modify the 



value of the fixed point in such a way that r* and u* become: 

This modification does not lead to variations of the eigenvalues of the matrix associated to 
the transformation, so that critical exponents are not changed from the white-noise case, as 
expected. 
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Hence we can conclude that the noises used here, either (p.6|) or ( |4.2| ), with a finite 
correlation time and length, change the position of the critical point but not the critical 
exponents, when compared to the white- noise case (|1.3| ). This is due to the way this kind of 
noises behave (see (p.3|), for instance). Should the correlation of the noise decay as a power 
law, then, according to the scaling analysis of Ref. [Q, critical exponents would be changed. 

V. COMPARISON WITH NUMERICAL RESULTS 

In Refs. model A ( |1 . 1|) with noises ( p..4|) and (|4.2|) was studied by means of a numer- 
ical simulation in 2-d. Two main results were obtained there. Firstly, it was established that 
either r or A stabilize the system, in the sense that the critical value of the noise intensity is 
enlarged. Secondly, critical exponents were also evaluated by means of a finite-size scaling 
analysis, and their values were found to be similar (within error bars) to the accepted values 
for the white-noise case. Since the numerical evidence of the critical- noise- intensity shift as 
a function of r and A is much clearer, that is what we want to explain in the light of the 
above DRG results. 

The simulation model, defined by fixed values of r, /i and u, was 

= ^ - V^^ + V^^) + e(x, t) (5.1) 

where the only independent parameters were the noise parameters D, r and A. 

Let us now consider Eq. ( p. . 1|) at the critical point with the parameters r*, u*, D*, fi* = 1, 
T and A. By means of a change of variables, one can recover the simulation model ( ^.1|) with 
a critical noise intensity given by 



D, = -^^ = De(r = 0,A = 0)^l + A^j +rj+0(r^) (5.2) 

where Dc{t = 0, A = 0) = = ^. 

Certainly one cannot expect a good agreement between a calculation in 4-d extended 
to 2-d and a numerical simulation of a 2-d model in a finite discrete lattice. In this sense. 
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the predicted value for Dc{t = 0, A = 0) is very different from tlie numerical result (0.38). 
Nevertheless, the ratio -Dc(t, ^)/ Dc{t = 0, A = 0) gives better results, as one can see in Fig. 
2. One can thus conclude that RG calculations have given a reasonable explanation of the 
critical results for model A when a colored noise is considered. 




Ratio of the critical noise intensity in the colored case to its value in the white 
case versus t (from Ref. ^). The lines correspond to the DRG result ^5.3{ ). 

Fig. 2 
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APPENDIX: CALCULATIONS OF THE DRG AT ONE-LOOP ORDER 



In order to eliminate the outer (high-momentum) modes from Eq. ( p.2| ) we break the 
integral term (the third diagram in Fig. 1) into modes located in the inner hypersphere (i.e. 
those with <| k |< Ae^'), which will be called external modes; and modes located in the 
outer shell (i.e. those with Ae"' <| k |< A), which will be referred to as internal modes. 



This leads to a rewriting of Eq. (|3.2| ), whose diagrammatic representation is shown in Fig. 
3 for both the internal and external modes. Fields and generators depending on internal 
modes are represented by means of a slashed line, and those depending on external modes 
are left unchanged. 



■ 

^ ^ — ^ — ^ 



^ ^ ^ 

^ ^ ^ ^ ^ 



(b) 

Evolution equations for the external (a) and internal (h) modes 
ready for the perturbation procedure 
Fig. 3 

The evolution equation for the internal field (Fig. 3b) can be solved iteratively up to 
(9(m^), and the result can be introduced in the internal contributions to the integral term 
of the corresponding dynamical equation for the external field (Fig. 3a). This leads to a 
renormalized equation where internal modes appear only in the noise term. After averaging 
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out this internal noise, the remaining renormahzed equation is the one represented in Fig. 
4. 
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Renormalized equation after internal noise averaging 

Fig. 4 

In order to understand the meaning of the new notation introduced in Fig. 4, let us 
make explicit the first diagram. It comes from an average over the internal noise term, as 
shown in Fig. 5. Explicitly, 

-uG<{k,uj) r r ^<(k2,cu2)(^°>(ki,a;i)^°>(k-ki-k2,^-u;i-^2)) = 

r> 2D 

-uG<{k,u;)^<{k,u;) ^^—G>{h,u,)G>{~k,,-u,) (Al) 

where definition of (|3.4| ) and correlation of the noise in Fourier space ( |3.6|) have been 
used. We will denote this result as shown in the second member of Fig. 4. The superindex 
< (>) on Go Ei'iid ip denotes an external (internal) mode. 



Interpretation of the diagram which renormalizes r 

Fig. 5 

Now we will be able to write down the transformation (recursion) relations for the pa- 
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rameters r, fi and u. We are in the zone where u is relevant {d < dc = 4) but are also 
supposing it to be small (which will presumably occur for ~ 4). Therefore if these re- 
cursion relations are to be an expansion around u = they should also be an expansion in 
e = 4 — d. Moreover, r will also be supposed to be small in the following calculations (we 
are interested in the critical point). 

It should be noted that the diagram in Fig. 5 contains only one thick line, and will 
therefore renormalize the diagram at the left member of the figure, whereas the second 
bubble contains three thick lines, so that it will add to the one-vertex diagram at the right 
member of the figure. That means, as we will see, that the first bubble renormalizes r and 
the second one u. On the other hand, it is easy to see that /x is not renormahzed at 0{u), 
since the one-thick-line bubble does not contain any term proportional to /c^. 

Concerning the computation of the diagrams, it should be said that the frequency in- 
tegrals can be evaluated by means of contour integrations in the complex plane, and the 
integrals in momentum space can also be calculated by assuming r ~ |^ and a very small 
T. These integrals should be expanded in e = A — d, but since they are all multiplied by u, 
which is also small, we only need their zero order in e. Thus we evaluate them in d = 4. 
We also let A = 1. Moreover, since the parameter / is the step of the renormalization pro- 
cedure, and we want this procedure to be continuous, we will assume / <^ 1. After all these 
considerations, the final result can be shown to be 

uD 

rj = r + 3 K^n- fir) I (A2a) 

uj = u-9 I A2b 

III = fi, (A2c) 

where is {27c)^'^ times the surface of the unit d- d-dimensional hypersphere: 

Kd = l^-'^Ti-'^l^ r(rf/2) (A3) 

We also want to know how the noise intensity D changes under the effects of the DRG. 
This is not given by the analysis explained above, because the parameter D does not appear 
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explicitly in Eq. ( p.2|) . To make it come out explicitly one can " autocorrelate" this equation, 
i.e. multiply the equation for ilj{k,uj) by the same equation for ilj{k',uj'). This can be done 
diagrammatically in a simple way. It can be seen that the first diagram renormalizing D is 
a two- vertex diagram, which means that it is of Therefore the recursion relation for 

D up to the order we are considering is simply 



D, = D 



(A4) 



It is worth noting that there is no first-order correction in r neither for /i, m nor D. 

Finally, we want the DRG procedure to lead to a renormalized equation as similar as 
possible to the original one. At this moment both equations are identical in form, the only 
difference being the range of the variable in momentum space to which they are applied. 
Indeed, the renormalized equation only governs the evolution of the external modes. In 
order to eliminate this difference, we shall perform the following rescalation: 



(A5) 



This scale transformation is identical to the one described in Sec. || with h = eK There 
we could see that the rescalation led to effective values of the parameters of the model 



and ( |2.4|) ). Hence, by combining the momentum-shell integration process and the 
space rescalation, one finds the following final discrete recursion relations for the model 
parameters up to first order in /: 



ri e ~ 



r + / 



uD 

z r + 3 (1 — /ir) 



(A6a) 



u = Uj 



^l(z-2a) r^^^l 



{z - 2a)u - 9 



(A6b) 



/i = /i/ e' 



(A6c) 



D = Di e'(2«+^-'^) c^D + /(2a + z-d)D 



(A6d) 
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